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ABSTRACT 



Quantized Frequency Modulation is one of the various 
methods used to reduce the effects of multipath propaga- 
tion of high frequency radio waves. The modulation 
technique is described, and a signal and noise analysis is 
made resulting in a theoretical maximum performance 
criteria. A comparison is made with the performance of 
various configurations of the system demodulator. 
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1. Introduction o 

It has long been recognized that reliability in high 
frequency radio communications is limited not so much by 
the power which can be radiated by an antenna as it is by 
the destructive vectoral addition of signals arriving at 
the receiver. The signals which are added vectorially 
may be delayed replicas or may be parts of two different 
symbols. In both cases the effect is the result of the 
delays in multipath propagation. The first case is selec- 
tive fading, and the second case is intersymbol interfer- 
ence . 

There have been many schemes devised in attempts to 
reduce the effects of multipath propagation phenomena. 
Several forms of diversity reception including time diver- 
sity, space diversity and frequency diversity each designed 
to isolate one of the incoming signals, and various 
combining techniques have been used, /l/ Other schemes 
include the use of single sideband /2/ and synchronous 
detection. /3 / 

Probably the most advanced technique in combatting 
multipath effects is the "Rake” system, /k / This system 
uses a cross correlation detection scheme, coupled with a 
method for isolating the signal from each path, and then 
combines the signals from each path in nearly optimum 
proportions . 

All these schemes have their own advantages and dis- 
advantages. Usually the best results are obtained by 
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those schemes requiring large antennas or complex equip- 
ment, or large land areas. Other schemes of less complex- 
ity can expect less increase in reliability. 

The modulation scheme analyzed in this paper is one 
designed to operate within the bandwidth of one single- 
sideband voice channel of 3 K.C. Also a minimum size and 
weight limitation is imposed which precludes large antennas 
or space diversity techniques. Instead a type of frequency 
diversity is employed which allows the transmitter to oper- 
ate at a nearly constant power output. A coded pulse is 
used for both the ’’mark" and "space", and a delay line 
summer is used to add coherently the output pulses. 

The modulation scheme will be described and related 
to the principles of communication theory. A signal and 
noise analysis will be made on the system to predict error 
rates for various signal to noise power ratios. 

The methods and results of laboratory testing of the 
demodulator under various signal to noise ratios will be 
presented. Comparisons will be made between the predicted 
error rates and the actual error rates. 
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2. System Operation. 

This section includes a description of the basic 
modulation scheme and the method used for demodulation. 

The principles of communication theory which this system 
implements are emphasized where employed. 

2.1 Basic modulation scheme. 

Quantized Frequency Modulation (QFM) is the name 
given to the modulation scheme developed by the Communi- 
cations Division of Hughes Aircraft Company. This scheme 
was devised in an effort to reduce the effects of multipath 
fading on the reliability of teletype or mark-space communi- 
cations systems. 

A combination of several techniques is applied, and 
these techniques will be described in later sections. 

The tapped delay line approximation to the matched 
filter is used in this system as follows. 

A combination of short pulses and frequency stepping 
is used to reduce the effects of intersymbol interference 
and selective fading respectively. The short pulses are 
spaced at each frequency to allow time for the path 
delays to die out. But upon detection these pulses are 
added coherently by means of the tapped delay line. Fig. 1 
is an idealized graphical representation of the signal 
transmitted. 
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Actual QFM Signal 
Figure 2 



This is the modulation on a carrier. Notice that there 
is uninterrupted modulation of the carrier, and hence the 
output power is constant. In order to contain the modulation 
within a 3&C bandwidth the pulses are shaped to approximately 
a cosine squared shape. An example of the actual shaped 
pulses is shown in Fig, 2 for Mark— Space-^Space—Mark as 
before , 

To demodulate this Quantized Frequency Modulation 
a demodulator of the type described below is used. Fig, 3 
is a block diagram which will be used in the following 
description. 

The QFM signal is obtained from a standard communica- 
tions receiver and is put through the four filters, but 
only one of four filters will pass a signal at any instant 
of time. The other three filters will pass only noise. 

The amplifier will amplify the signal to a maximum level 
if the signal is large enough to reach that level, other- 
wise it is amplified by the full gain of the amplifier. 

Each channel is full wave detected either positively or 
negatively, and the cnannels are added as shown in Fig. 3» 

The signals are then summed in the tapped delay line and 
differential amplifier. 

The output of the differential amplifier is used to 
establish synchronization, which in turn gates the inte- 
grator to operate only during the time that the peak signal 
occurs . 
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The output of the gated integrator will be positive 
or negative depending on whether a ’’mark” or ’’space" code 
was transmitted o The output of the gated integrator is 
sampled just before the integrator capacitor is discharged „ 
The sampling timing is also controlled by the synchroni- 
zation circuit o The polarity of the integrator output will 
then determine v/hether the data is ’’mark" or "space" and 
\\rill trigger the appropriate data flip-flop „ 

Since the decision is based solely on the polarity of 
integrator output the decision is made about a nearly 
zero threshold o 

The operation of the system has been described „ The 
principles of communication theory employed will be described 
where applicable in the following sections 0 
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2,2 Short pulse and multi-frequency transmission 

The effect of multipath delay for pulses which are 
long compared to the delay time is shown in Fig 0 4, 

The effect of multipath delay for pulses which are 
short compared to the delay time is shown in Fig 0 5* 

It can be seen that the short pulse is more desirable, 
but of course this pulse requires more bandwidth, and con- 
tains less total energy 0 

Turning now to the selective fading effects of multi- 
path propagation and directing our attention to Fig, 6, 
it may be noted that the multipath delay is different for 
different frequencies. So one possible approach to com- 
batting selective fading could be by using more than one 
frequency in the transmission. 

In this QFM scheme, the techniques and those to be 
described in Sections 2,3 and 2,4 are combined. In 
addition, in order to increase the total energy in a mark 
baud and space baud, ten pulses are transmitted for each 
baud. At the receiver the envelopes of these pulses are 
added coherently as described in the previous section. 
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2«3 Coded pulses and matched filter detection 



It has been shown / 5/ that the optimum receiver to 
detect a signal in white Gaussian noise is one which 
calculates the likelihood ratio and then decides if a 
signal was present or not depending on whether the like- 
lihood ratio exceeded a given threshold value or not 
respectively,, The receiver which does this is a cross 
correlator and essentially calculates the coefficient of 
linear correlation between the received signal and a copy 
of the transmitted signal 0 If this exceeds a specified 
value then the decision is made that a signal was present „ 
Since this is a continuous process C(T)is continuously being 
calculated as 



where is the received signal and 

S^ft) is the copy of the transmitted signal 0 
For application to the communications problem the 
decision to be made is, "Was the signal sent a mark or 
a space?" c It has been shown, /l/ that the optimum 
decision process is one which calculates for both 
the mark and space by cross correlating the incoming 
signal with the "mark" signal as transmitted and the 
"space" signal as transmitted, and then decides which 
is larger,, 




(o 



Thus 
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and ^nyut'fc- (?) ' tPoo ~=F 



l r/ ^ (t-r) At & 

where Ssf«tt {O is the signal representing "Space" 
and Snu>rM.(t> is the signal representing "Mark". 

It has also been shown /l/ that the probability of 
making a wrong decision is least when mark space 
and the decision threshold is zero. The cross correlation 
of the "mark" signal as transmitted with a "space" signal 
as transmitted for this case would be -1. 

One means of implementing a cross correlation detec- 
tor is by means of a filter. The output of a filter for 
a given input signal is the convolution of the signal 

with the impulse response of the filter. 

r o* 

' J X (?) k 7) A T (V 

Or 

P* ' J-a *(t> k fc-t) At (f) 

if the impulse response of the filter h(t) is the time 
reverse of the input signal, X(+t), then the output is the 
convolution of X(t) with X(-t). Substituting in equation 
( 5 ) above 

CO 

# (Tj *(£> * (T- I- *i) At (i) 

£ flv '-J* x%) x (jt t) At C 7 > 

P CO z x ft* V ^ ( e > 

-<x> 
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Since \(y~) as shovm in equation (1) is an even function 
for periodic signals 

fnv = C(ci) (V 

■f r c yh (i ) 

(V; ((-t> - r^> ' / \o * if. v & (,c> 

J - T A 

In any physical implementation, integration starts 
at some time t 0 and is carried out for some finite time 
interval. Hence 

Ctr> 

for any physically realizable system. 

In terms of probability, the question, "Given that 
a signal was sent, what is the probility that it was a 
"mark” and what is the probability that it was a "space” ?" 
must be answered. Since each is equally likely, the a 
priori probability of each is one-half. 

The a posteriori probability that a mark (for instance) 
was sent is, from /l/ 

? -• - r 

f/fp/s +*i*rfa {t)~] flS^UrK (tj~] +?[{(£)/■&/’" •-‘-(t'] P (t) J 



It has been shown /l/ that for mark and space equally likely 
the decision reduces to: 
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"Decide "mark" if 

r r s T / 

J 0 fit) ark (t) <rf J Q f(t)Syi/c.t(£)df' h Z ) 

is greater than zero (positive). Decide "space" if (12 ) 
above is less than zero (negative). 

If a truly matched filter or correlation detector 
could be built the error rate has been shown to be 





Figure 7 



2.4 Tapped delay line approximization to a matched filter. 

Turin /6/ proposed three methods of implementing a 
matched filter. One of these was the tapped delay line. 
However, he pointed out that there are physical limitations 
to this approach. The filter which is truly matched is 
matched continuously from beginning to end of the signal 
when T = 0. Fig. 8 shows an example of a signal and the 
impulse response of a filter matched to it. 




signal 



h(t) for filter matched to signal 
Figure 8 . 

A signal which might be put into a tapped delay line and 
its corresponding matched filter is shown in Fig. 9. 
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Figure 9 

This particular signal pattern chosen is one of the 
well known "Barker codes" /7,$/. It was chosen to illus- 
trate a desirable property which is possessed by some 
binary codes, namely that of a "good" autocorrelation 
function. Notice that the interval between peaks is flat. 

However if instead of a truly matched filter a 
tapped delay line is used, the output of the filter is 
still the correlation of the input signal with the impulse 
response of the filter, but the taps are not as "broad" 
as the pulse length. Instead the taps could be represented 
as delta functions weighted by the proper sign as shown in 
Fig. 10. 
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Figure 10 

Notice that the output^ at T = 0 has the shape of 
an individual signal of width T instead of the triangle. 
The delay line it is seen, linearly adds the voltages at 
each tap, and is not truly "matched” to the signal. The 
next section will deal with the analysis of the system, 
and will establish a mathematical model to predict theoret 
ical performance maxima* 
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3. Analysis of Demodulator. 

3.1 Mathematical model of demodulator. 

In order to analyze the performance of the system 
and to be able to predict error rates for various signal 
to noise power ratio inputs a mathematical model was 
devised. Models have been established and performance 
analysis done for several of the well known binary 
communications systems. /9,10,ll/ However most of these 
analyses assume a signal to noise ratio of 3db or more. 

This paper deals with lower signal to noise power ratios 
and consequently the approximation used in them are not 
valid for this study. The signal to noise power ratios 
in and out of an envelope detector have been described 
completely by Rice. /l2/ Techniques in the application 
of his results may be found in other publications as well. 

/13, 14/ 

The exact expression for the density function of the 
output of an envelope detector includes a modified Bessel 
function. At this point the high speed digital computer 
was used to obtain the signal voltage and noise power for 
different signal to noise ratios. This system is essentially 
linear from input to the full wave rectifier of each channel, 
and subsequent to that. Hence if the non-linearity of the 
full wave rectifier could be analyzed and modeled then the 
entire system could be modeled. The amplifier prior to 
detection will be assumed linear, and section two demon- 
strated that the tapped delay output was the linear alge- 



braic sum of the signal in the filter at every instant 
of time. Hence the signal to noise power ratio at the 
input to the full wave detector will be taken as a refer- 
ence. Signal to noise power ratios at the input to the 
demodulator can then be expressed in terms of the refer- 
ence. In making the transition from full wave rectifier 
input to demodulator input it will be assumed that the 
noise figure for the AGC amplifier is zero. 

Another assumption in the model is that the signal 
pulse is a square pulse. This is done in order to simplify 
the analysis of the model. However in order to use the 
model to relate predicted error rates to actual error 
rates measured a conversion factor is used. Appendix IV 
shows in detail how this conversion factor is calculated. 
Once the relationship is made the relative positions of 
the predicted curve to the actual measured error rate curve 
is established. The error rate versus signal to noise 
ratio per pulse at the input to the detector can be used 
as a reference, and error rates for signal to noise ratio 
at any point before the detector can then be obtained. This 
relationship may be seen by comparing the graphs of Fig. 

19 and Fig. 20. 

The model is considered to have four separate channels. 
Each channel has an ideal filter 580 cycles wide centered 
at 680, 1260, 1940, and 2620 cps. The signals of each 
channel are envelope detected and added together as shown 
in Fig. 3 • The output of the tapped delay line is inte- 
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grated for one half the pulse period during the time 
interval that the summed pulse is highest. A decision is 
then made at the end of integration time on the basis of the 
polarity of the integrator output. "Mark" if the integra- 
tor output is positive and "Space" if the integrator output 
is negative. 

The probability of error is the probability that 
the integrator output will be negative (or positive) when 
a "Mark" (or "Space") signal was transmitted. 



From Fig. 11 it may be seen that the probability of error 



In order to take advantage of the tabulations of the 
area under the normal curve /15/ this is put into the form 




Figure 11 



is : 







where 
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3.2 Analysis of model. 

It has been shown by Rice /ll/ that if an envelope 
detector has as its input bandlimited white Gaussian noise, 
the probability density function for the output is the 
Rayleigh distribution. Specifically for an input or 
bandlimited white Gaussian noise of power spectral density 
N 0 watts per cycle per second, the output of an envelope 
detector has as its probability density function the Ray- 
leigh distribution, 

fs (y)= 1_ ye-y 2 /2 N 2 

f N 2 

where N 2 = N 0 B. 

The power spectral density and probability density 
function is as shown in Fig. 12. As a signal is added to 
the noise input, the density function shifts to the right 
as shown in Fig. 13 and Fig. 14. 

The probability density function for the output is 
now a function of the signal power and the noise power. 
Specifically from /l4/ 



where I Q (x) is the Modified Bessel function of the first 



signal to noise power ratio for a continuous carrier. 
Most treatments of envelope detectors will make a closed 




kind and zero order for the argument x, and s 2 is the 
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form approximation to I 0 (x) for large or small signal to 

noise ratios, because the density function becomes 

essentially normal for large signal to noise ratios, and 

2 

I (x) can be closely approximated by e£ for x<<1. How- 

4 

ever for this study these approximations were crudest in 
the region of main interest. The modified Bessel function 
I Q (x) may be expressed as the sum of an infinite series, 
namely 



oO 



<x> « Z 



Zh 






Appendix I shows the infinite series solution to ob- 
taining the mean, and variance of the resulting probability 
density function. 

Appendix II is the computer programs used to calcu- 
late the mean and variance of the distribution resulting 
from different signal to noise power ratios at the input 
to the envelope detector. The program also includes some 
other calculations used to obtain the probability of error 
which would have been quite laborious if done by hand. 

Also included in Appendix II is a sample of the program 
outputs . 

Fig. 16 illustrates an example of the probability 
density functions for all four channels at one instant 
of time, and how they are combined. 

Although the distribution of the sum of the four 
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channels is not truly normal, it is nearly so. After 
summing ten of these nearly normal distributions, by the 
central limit theorem, the resulting distribution will be 
very close indeed to being a normal distribution. 

This signal is then integrated for an interval of 2 
milliseconds. Appendix III describes the treatment of 
integrating a signal plus noise. For a continuous carrier 
pulse this amounts to reducing the signal to noise power 
ratio. The result is manifest on the distribution function 
as a change in the ratio of the D. C. signal level to the 
standard deviation of the distribution, and the probability 
of error is reduced accordingly. Probability of error is 
now 



Pg * 2 1 -Erf Q 

/ o {cM i ~ 1- ) (/<■?) 

where Q = — — 

JL ( is the mean of distribution of the output voltage of 
the channel containing the signal during one pulse dura- 
tion. is the variance of the distribution of the 

output of the channel containing the signal during one 
pulse duration, or fluctuation noise poi^er. 

is the mean of the Rayleigh distribution resulting 

_ 2 

from envelope detection of white Gaussian noise, o 2 ^he 
variance of the Rayleigh distribution or fluctuation noise 
power of the output of envelope detection of white Gaussian 
noise. 
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PROBABILITY OENSITY /»(v| 




Probability Density Function for the Output 
of an Envelope Detector for Increasing 
Signal to Noise Ratio. From /II/. 



Figure 1 5 



v 
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T 



The standard deviation of the Gaussian noise distri~ 
bution was taken as one, and thus the noise power was 

one. Then the signal to noise power ratio is just the 
signal power. A Control Data Corporation 1604 computer 
was used to calculate Q and the probability of error 
obtained from tables /15/. The probability of error for 
various signal to noise power ratios was obtained and 
plotted. Curve V of Fig. 19 is the curve obtained. 

Laboratory measurements were made to compare actual 
error rates with those predicted. A description of the 
apparatus and techniques used and the results obtained is 
contained in the next two sections. 



* 
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4. Laboratory Tests of Demodulator, 

A laboratory test of the demodulator was desired to 
compare actual performance with the theoretical predictions. 
In order to do this it was necessary to devise some means 
of detecting and counting errors and counting decisions 
made. Appendix V is a description of the error detector 
which was built up from Engineered Electronics Company 
digital building blocks. Fig. 17 is a block diagram of the 
laboratory equipment layout. 

The automatic gain control circuits of the amplifier 
were set for minimum gain, and the linear range of the 
amplifier was determined. Fig. IS is an example of the 
linear range of the amplifier. 

A few tests showed that the error rate was quite 
good for good signal to noise ratios but seemed to have 
a signal to noise ratio threshold below which the demodu- 
lator lost synchronization. Actual errors as a result 
of wrong decisions were still quite low. 

Investigation of the synchronization circuits re- 
vealed that for the signal to noise ratio at which the 
demodulator lost synchronization, there was still available 
a reasonably "clean" signal pulse coming from the differ- 
ential amplifier. However the voltage level of the signal 
was not of sufficient amplitude to be sensed by the part 
of the synchronization circuit which decided if a signal 
(either "mark" or "space") was present. A two stage ampli~ 
fier was then inserted between the differential amplifier 
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and the "signal sensing" portion of the synchronization 
circuit. Curve II of Fig. 19 is the curve of error rates 
measured for this configuration. 

Visicorder recordings of the outputs of the gated 
integrator and data output showed that additional improve^ 
ment could be obtained. The time constant of the operational 
amplifier integrator was reduced. The error rates obtained 
for this configuration is shown as curve III of Fig. 19 ° 

Visicorder recordings show that the sensing threshold 
was considerably above zero. The input resistance to the 
sensing transistors was reduced to allow the required 
switching current to be reached at a lower threshold 
voltage. This was the final configuration tested, and the 
error rates for it are shown as curve IV of Fig. 19° 

In Fig. 19 the measured signal to noise ratio was 
corrected to fit the model as described in Appendix IV. 

Fig. 19 is the same data referring the model signal to 
noise ratios to the measured signal to noise ratio per 
bit. 

As discussed in section 3, the data in Fig. 19 and 20 
are for the signal to noise ratio at the input to the full 
wave rectifier as a reference, in order to refer the error 
rate to the signal to noise ratio oer bit at the input to 
the filters the following correction is required. 
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Each filter bandwidth is 5^0 cps. The total noise 
bandwidth at the input is 3 kc — 400 cps = 2500 cps. There 
is then or 4.5 times as much noise per bit at the 

input as in each channel. Signal to noise ratio must then 
be reduced by 6.53db. Fig. 21 shows this. 

Any measurements made in the future on this system 
will probably be done by a true reading RMS voltmeter as 
was done in our experiments. Figs. 19, 20, and 21 are on 
a signal to noise power per bit basis. Since the signal 
was a mark — space — mark — space, etc., the average signal 
power was just one-half of the power per bit. Average 
signal to noise power ratio is then 3db less for each 

rate than for Fig. 21. This is shown in Fig, 22, 
and is the error rate for the signal to noise ratio which 
would be measured at the input to the demodulator. 

Fig. 23 is a group of pictures showing the QFM signal, 
the differential amplifier, output for two signal to noise 
ratios and configurations. 
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Probability of Error 
vs . i 

Average Signal to 
Nol3e Power Ratio 



at Input to D emodu 



lator, Referred to a 
Cosine Squared Pulse 



36 





5v/cm .J'K signal (no noise) 
5v/cn. differential amplifier 
output 




lOv/cm I»i signal + noise 
3/u (per bit) = 243db (Ref COS 2 ) 
Iv/cm output of differential 
amplifier 
P E (II) = .16 




5v/cm gated integrator output 
Iv/cm output of differential 
amplifier 
P E = . 16 ( 11 ) 




,5v/cm differential amplifier 
output 

S/K (per bit) = .3db (Ref COS 2 ) 
5v/cm gated integrator output 
PE (IV) = ,o 




5v/cm gated integrator output 
5v/cm data outout 

P E (IV) =^0*7 



F 3 3 









5. Conclusions. 

A close analysis of a system to obtain a theoretical 
optimum performance level provides a useful standard to 
which actual performance may be compared. Such an anal- 
ysis may show the possible improvement to be gained by 
alterations in various sections of a system. Such was 
the case for the QFM system studied. Curve V is the 
theoretical optimum performance criteria for this system* 
and I to IV are performance levels for various configuations 
of the system. It appears that additional improvements in 
system performance are still possible. However these 
curves cannot show the limitations of the synchronization 
circuit. Perfect synchronization and zero threshold are 
assumed. 

In any physically realizable system synchronization 
is not perfect and the threshold is some small value above 
and below zero. 
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APPENDIX I 



Infinite Series Solution to Mean and Variance of the 
Probability Density Function for the Output of an 
Envelope Detector for Various Signal to Noise 

Power Ratios 



From /ll/ the probability density function for the 



output of an envelope detector is 



-S' 



- Al 



fp-)-- c ^ x o(ks^L-) 

where 

S^ is the signal to noise power ratio at the input 
to the envelope detector $ 

N 2 is the noise power at the input and 
I 0 (Y) is the modified Bessel function of the first kind 
and zeroth order of the argument Y « 

I q (Y) can be expressed as the sum of an infinite series 

I. (Y) Z „ 

1 * (//■'} 

The mean of this distribution is 
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since the distribution has only positive values for X and 
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Term by term integration yields 
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For digital Computer Operation this series may be 
nested as 








r" ~ — ■ " — ”p — — T 

* ^ -? • 3 V ;> • y v 



The computer program which does this is in Appendix II „ 
The Second Moment for the distribution function has 
been evaluated in closed form as 



b [ x*-] = * 

The variance or fluctuation noise power is then: 

cT^ - EWI -{£l xlT 

^ ' R 2 an d ^may be evaluated for different signal to noise 
power ratios. 
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Appendix II 

The digital computer program written to compute the 
mean and variance of the distribution for the output of 
envelope detector is shown below. Sample outputs are also 
included in this appendix. 



PROGRAM SIGNOIS 

SIGNAL VOLTAGE AND NOISE POWER FROM AN ENVELOPE DETECTOR 
, D ££PV5?8!i. T 5 M U o:> >A EM1 (90),EM2(90) ,X0ISPR(90> ,SIGPR<90) 

1 ,DCSIG(90) ,EM1SQ3(90) ,TDCSIG(90),TN0S PR( 90) ,ERF ( 93 ) 



• VJ / ten low J i • i iu ' 

2,DB(90) ,SIGINT(90),ERFINT<90) 
P0WER=0 
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DO 20 1=1,90 
1 I F { 1-10)4,4,5 

4 P0WER=P0WER+.01 
GO TO 10 

5 IF ( 1-40)6,7,7 

6 POWER “POWER*. 1 
GO TO 10 

7 IF( 1-60)8,9,9 

8 P0WER=P0WER+.5 
GO TO 10 

9 P0WER=P0WER+1 .0 

10 SIGPR( I )=P0WER 
XEGPR=-SIGPR( I) 

T= 1 • 0 

SUM= 1 . 0 
TRM ( 1 )=1.0 
DO 15 N=2 ,100 

T=T+ 1 . 0 , 

0TRM(N) = ( (TRM(N-l) ) * ( S I GPR (I) * ( 2. 0*T-1 . 0 ) ) )/(2.0*( I T-1 .0)*/(T-l .0) ) 
I SUM=SUM+TRM(N) ' 

EM1(I)=SUM*EXPF(XEGPR)*1. 2533141 
EMI SQD( I ) =EM1 ( I )* EMl ( I ) 

EM 2(1 )=2.0*(1.0+SIGPR(I) ) 

XOISPRl I ) =EM2 ( I )-EM1 SQD( I ) 

DCSIG ( I ) =EM1 ( I )-l .2533141 
TDCSIG( I)=10.0*0:SIG( I) 

TN0SPR( I ) =10.0* (XOISPRI I ) *1.28761 102) 

QPR=TNOSPR( I) n 

20 ERF( I )=TDCSIG( I )/ SQRTF (QPR)l 
- PRINT 21 

21 OFORMAT ’(97H SIG TO NOISE IN FIRST MOMENT OUT 

1 ? OUT DC SIGNAL VOLTAGE NOISE POWER OUT) -> 

PRINT 3, (SIGPR(I) , EMI (I) , EM2 ( I), DCSIG ( I ) , XOISPR ( I ) , 

1 1=1,90) 

3 FORMAT (5F20.8) 

PRINT 29 

PRINT 98, (SIGPR(I ) , TDCSIGt I ) ,TN0SPR( I ) , ERF( I) ,1=1,90) 

98 FORMAT (4F20.8) 

PRINT 29 

29 FORMAT! 1 H 1 ) 

DO 30 1=1,90 

DB(I)=10.0* L0G13F ( SI GPR ( I ) ) 

SIGINT(I)=1.2*TDCSIG( I) 

30 ERFINT ( I )* S IGINT ( I ) /SQRTF ( TNOSPRC I ) ) 

PRINT 31, (DB(I),SIGINT(I ), ERFINT ( I), 1*1, 90) 

31 FORMAT (3F20.8) 

END 

END . , 

PRINT 99, (TRM( I), 1-1, 100) 

99 FORMAT (F20.8) 



SECOND MOMEN 
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APPENDIX III 



FILTER CHARACTERISTICS OF AN INTEGRATOR 
The impulse response of an ideal integrator which 
integrates for a time interval T is shown as a solid line 
in Fig, A, /17 / 
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Impulse response of an ideal integrator. 
Figure III-A 



To obtain the transfer characteristics in the frequency 
demain translate the coordinates an amount T/2, shown as 
dotted in Fig, A. Transform this to the frequency domain. 
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A delay in the time domain of T/2 seconds is manifest 

-4 i*J TA 

as a phase shift in the frequency domain of <£ ^ 0 Hence 

H(^> ~ <f ^ T/i_ 



= ^ rr £r_ 

rrf-r 



The power output of such an integrator would be 
r cc ^ 

/$ - / ) I 6-(lj) 4 -f 

^-00 

where G(W) is the power spectral density of the input 
signal 0 

The power output of the integrator becomes 

■CO 

r n = / I 
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for one sided frequency spectrum 0 

For this problem T ~ 2 x 10 seconds „ The frequency 
of the first zero is then 

rr { r -- tt 

■f = yT “ S~ 0 0 C,f> ,S . 
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For any one of those three channels which for any 
instant contain noise only the power spectral density is 
as shown in Fig. 13 from /ll/ 
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From the Rayleigh distribution, the variance or fluctuation 
noise power is ,429tfN 2 . 

Then — 
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There will always be three channels of noise power 
,4290'fj 2 and the noise power from the signal channel going 
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into the integrator which will add directly since the 
noise from different parts of the spectrum are not coherent 0 
The noise power in the signal channel will be increasing 
with the increasing signal power e Thus the noise power 
may be considered as four channels of noise power plus 
the increased noise power in the signs 1 channel due to 
the signal times noise products. This increase may be 
closely approximated by curve fitting techniques applied 
to Fig, III-B. By linear regression the increase in 
noise power is 



where X is the signal to noise power ratio. 

This increase in noise power is essentially flat over 
half the original bandwidth of the noise input to the 
detector. Thus the Noise Power Spectral Density into the 
integrator will be 




G-(w) - 

^ ri n 



s ~0 O 




and the power output of the integrator is 




.too 
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The power spectral density input to the integrator- 



will be as shown in Fig. III-C. 




In order to obtain a tractable model the noise power 
spectral density will be assumed to be flat from zero to 
B/ 2 cycles per second with a total noise power equal to the 
original noise power. The noise power out of the detector 
then is 




Transforming variables 

r /T /- 

Fo - ^ J & - . ? ? ^ y 

* & 



This modifies the error rate to be 

l o ( , ~ ^ 

</-;(' ^(<*7 v 




i r __ /- / 0 ( ^ v ) I' ^ ~1 

-r l-Erf — ■ - — 

z *- / / O ( ^7 *> 3 V J 
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APPENDIX IV 



The analysis of the QFM demodulator assumed a 
square pulse shape for each bit transmitted as shown in 
Fig. 1. The pulse transmitted is shown in Fig. 2. In 
order to bridge the gap between the model and the actual 
modulator the following conversion proeeedure was applied. 

The power in a constant unmodulated carrier of 
amplitude A is A 2 /2. 




4 



n 



Pc = 



2 , 




The cosine squared pulse which is the output of the 
rectifier is obtained from an input to the rectifier of 
a waveform like a carrier 100 per cent modulated by a 



sine wave. The power in this signal is 3/2 



A 2 



In 



terms of the peak voltage of the modulated carrier 






3 [k. - 2 - Cc 

2 Sr t 2 
The ’’mark" or "space" decision is based on the polarity 



2 (r^ ° 



of the integrator output after 2 millisiconds of integration. 



The output of the integrator for the modulated input is 

'tz / /* £ o ' 



T A 



% 



The amplitude of a constant carrier pulse for an equivalent 
voltage output of the intigrator would be .82 PK. 
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The power in this carrier is 

£? _ (• Z A Pk)~ ,6?2V- Fa 

r c. r - ■ r — 



* 

The power in a constant carrier signal which will produce 
the same output at the integrator as a modulated carrier 
with a power per bit of Pg is 



F E 4 r P & 



• z. Z 3 A 



Vs 

1.7 Z Fa 






Since all signal to noise ratios are referred to a noise 
power of one, the Signal Power is the signal to noise 
ratio, and 

ID Lop Pp ~ S 0 L op l ?? 1 10 
^//y £(j " q ^ 



Thus 2c53db should be added to the signal power per bit 
as measured. 

To refer the model to the measured signal power 2„53db 
must be subtracted from the theoretical noise power. 
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APPENDIX V 



ERROR DETECTION INSTRUMENTATION 
Some means of detecting errors was necessary in 
order to obtain an error count . Since the output data was 
digital, a digital system was used for error detection., 
Engineered Electronics building blocks were arranged as 
shown in the diagram below. 




Flip flop A assumed the state corresponding to the "Mark” 
or "Space" data received. Flip flop B assumed the state of 
the transmitted data. The outputs were gated according to 
the following equation. 

C * (A . B) - (A . B) 

Error = ( C , D ) 

D was derived from the interface timing circuits. 
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